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I .  INTRODUCTION 

In  1955,  four  shell  types  exhibited  unusual  behavior  which  involved 
the  movement  of  their  internal  parts.  Three  of  these  shell  - the 
20mm  M282E1,  the  20mm  T216E1  and  the  3Cmm  T306E10  - showed  much  less 
fast  mode  damping  when  fired  with  the  M505  fuze1"3.  This  fuze  has  a 
spherical  arming  rotor  in  a cylindrical  cavity  with  small  but  non- zero 
clearances  (Figure  1).  The  fourth  suell  - the  8-inch  T317  - showed 
significant  range  losses  and  very  large  spin  decays4.  This  shell  has 
several  rings  held  on  a central  column  but  free  to  move  with  small  but 
non-zero  clearances.  In  all  these  cases,  small  amplitude  motions  of 
the  internal  parts  had  significant  effect  on  the  parent  s> ell's  motion. 

Two  types  of  internal  ^arts  motion  are  possible:  (1)  linear 

movements*  of  their  centers  of  mass  relative  to  the  external  shell 
center  of  mass;  (2)  angular  motion  of  their  spin  axes  with  respect  to 
the  external  shell  spin  axis.  In  order  that  these  internal  movements 
have  a significant  effect  on  the  shell  motion,  they  must  have  components 
at  the  same  frequencies  as  the  shell's  pitching  and  yawing  motion. 

Thus  we  have  an  internal  resonance  situation,  but  the  amplitude  of  the 
internal  movement  is  bounded.  In  this  report,  we  will  develop  the 
theory  for  internal  motions  that  have  fixed  phase  with  respect  to  the 
angle-of-attack  plane  and  obtain  quasi- linear  solutions  to  the  resulting 
differential  equations.  Next,  much  simpler  results  will  be  obtained  for 
motion  which  has  fixed  phase  with  respect  to  the  plane  of  the  higher 
frequency  component  of  the  pitching  and  yawing  motion.  Finally,  these 
results  will  be  used  to  explain  the  observed  behavior  of  the  four  1955 
shell . 


*Uodapps  has  derived  general  equations  for  this  motion  but  applied  them 

to  the  very  simple  case  of  longitudinal  motion  of  a mass  on  the  shell ' s 

axis  of  symmetry . 

1.  E,D.  Boyer , "Comparison  of  Aerodynamic  Characteristics  of  20mm  HEI 
Shell  M97  with  Fuze  M75  and  20mm  Shell  T216E1  with  Fuze  M505," 
Ballistic  Research  Laboratories  Memorandum  Report  865 , \pril  1955. 
AD69009. 

2.  E.D.  Boyer , "Aerodynamic  Characteristics  for  Small  Yaws  of  20rm  Shell , 

HEI , T282E1  with  Fuze  M505  for  Mach  Numbers  . 36  to  3.  78, " Ballistic 
Research  Laboratories  Memorandum  Report  916 , August  1955.  AD77515. 

3.  E.T.  Roecker  and  E.D.  Boyer , "Aerodynamic  Characteristics  of  30mm  HEI 

Shell , T306E10,"  Ballistic  Research  Laboratories  Memorandum  Report 
1098 , August  1957.  AD152952. 

4.  B.G.  Karpov  and  J.W.  Bradley , "A  Study  of  Causes  of  ":ort  Ranges  of 
the  8"  T317  Shell,"  Ballistic  Research  Laboratories  Report  1049, 

May  1958.  AD377548. 

5.  A.E.  Hodapp,  "Equations  of  Motion  for  Constant  Mass  Entry  Vehicles 
with  Time  Varying  center  of  Mass  Position Sandia  Laboratories 
SC-RR-7 0-691,  November  1970. 
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II.  THEORY 


The  theoretical  model  consists  of  two  bodies:  (1)  the  external 

symmetric  shell  body  with  mass  and  (2)  the  internal  rotationally 

symmetric  component  with  mass  m . The  center  of  mass  of  the  internal 

component  is  allowed  to  perform  a circular  motion  of  radius  e normal 
to  the  axis  of  symmetry  of  the  external  body  maintaining  a fixed  phase 
angle,  <j>£,  with  the  angle  of  attack  plane.  The  axis  of  symmetry  of 

the  internal  component  can  cant  at  a small  fixed  angle,  y,  with  respect 
to  the  axis  of  symmetry  of  the  external  body.  The  plane  of  this  cant 
angle  is  now  constrained  to  rotate  about  the  shell  axis,  maintaining 
a constant  angle  $ with  the  angle-of-attack  plane.  If  0 is  the  roll 

orientation  angle  of  the  angle-of-attack  plane  with  respect  to  aero- 
ballistic  axes  , the  polar  angle  of  the  circular  center  of  mass  motion  : 


6 

E 


= 0 t d) 

Te 


and  the  roll  angle  of  the  cant  plane  is 


(] 


0 

Y 


0 + d> 

Y 


(2 


In  Appendix  A,  the  angular  momentum  of  this  two-body  system  is 
computed  and  used  to  obtain  the  following  differential  equations. 


*xbPb  + *xcPc  AJt  Pb 


- B (3  cos  0 + a sin  0 ) 

Y Y Y 


.[« 

P 


av 


B (6  cos  6 + a sin  0)0  e 


0 


av 


(A  + iLx0)C  -(Aa  ♦ iA  K + i (f  ♦ E)  = 0 


pa' 


(3 

(4 


6.  C.H.  Murphy,  " Free  Flight  Motion  of  Symmetric  Missiles,"  Ballistic 
Research  Laboratories  Report  1216,  July  1963.  AD442757. 
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where 


\ = 6 + i a is  the  complex  angle  of  attack 


I , , 1 are  the  roll  moments  of  inertia  of  masses  b and  c 
xb  xc 


Itb>  Itc  are  the  pitch  moments  of  inertia  of  masses  b and  c 


B = I p - I.  6 
Y xc  rc  tc 


B = m x l 
e c c 


r = B^  y exp  (i6  ) 


E = B£  0 e exp  (i0£) 


^xO  = *xb  ^b  + *xc  Pc 


and  the  other  symbols  are  defined  in  the  List  of  Symbols. 

The  usual  quasi- linear  assumption  is  made  that  the  pitching  and 
yawing  motion  when  the  internal  motion  occurs  has  the  same  epicyclic 
form  that  describes  the  rigid  body  motion: 


i$i  id>2 

5 - Kj  e 1 + K2  e 2 


where 


K.  = A.  k. 
1 J 1 


= <J>.  + $.  t 

In  Appendix  B,  quasi-linear  relations  for  the  frequency  shifts 
and  damping  effects  introduced  by  these  internal  motions  are  given. 


4 V - *jr  “J  Cj  I2  :t  +jr  - Lxo) 

K.  . XjK.  * S.  I21t  - Lx0]“ 
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where 


c. 

J 

- [B  .y 

L YJI 

cos  d>  + 
Y 

s. 

II 

t-J. 

s in  4>  + 

J 

Y 

B . 

= I p 

- I + . 
tcj 

YJ 

xcrc 

B $.e  cos  <{>  ] [e*9] . 
eTj  TeJ  L 

B | .e  sin  4>  ] [e10] . 
erj  eJ 1 


d4> 


4>  = 4>  - $ 

2 l 

xj  * [Vi  * V[2V>  • W 


and  the 


are  the  frequencies  of  the  rigid  body. 


The  quasi- linear  approximations  can  be  inserted  in  the  roll  equation 
(Equation  3)  and  the  result  averaged.  This  is  done  in  Appendix  C 
with  the  result 


IuPi+I  P ■A,  a - HS  - <j>  K S 
xb^  xc*c  l T>  Ti  i i y2  2 2 


(8) 


It  should  be  noted6  that 


<t>  + 4 « L J I 

Tjr  T2r  xo'  t 


(9) 


Thus  21  | . 

^ J 

slow  rate. 


L is  positive  for  the  fast  rate  and  negative  for  the 
xO 


Equations  (7-8)  show  that  positive  S.'s  ($  ,6  in  the  first  or 

3 c Y 

second  quadrant)  can  cause  undamping  of  the  fast  component  (K  ) , damping 
of  the  slow  component  (K  ) and  a reduction  in  the  spin  rate. 
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saipfps 


Now  0 can  be  determined  by  the  definition 


where 


Thus 


where7 


6 = 


r i0i 

[e  ] 


[6"1] 


e2 


i6  ~ K-\ 
e = 5 6 1 

(10) 

ki  = (K  2 + K 2 + 2K  K cos  4>) 1//2 
1 2 12 

K ^ 

2^  vf  [1  + 4 

„ 2tt 

h A /v  /v 

^ [1  + (K  /K  ) cos  <t>]5_1  d(j) 

w 2 1 

0 

K [6-1] 
ll 

(ID 

K Kl]e, 

2 e2 

(12) 

— — • [ (K  -K  ) E (k)  + (K  +K  ) E (k)]  (13) 

rrK  2 12  1 12? 

1 


— [(K  -K  ) E (k)  + (K  +K  ) E (k)]  (14) 

uK  2 2 11  12?. 

2 


4K  K 

k2  = L_2 — 

(K  +K  )2 
1 2 


7.  C.H.  Murphy,  "Prediction  of  the  Motion  of  Missiles  Acted  on  by 
Non-linear  Forces  and  Moments Ballistic  Research  Laboratories 
Report  995,  October  1956.  AD122221. 
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and 


E^  is  tne  complete  elliptic  integral  of  the  i-th  kind. 
When  K^>K^,  Equations  (13-14)  can  be  approximated  by 


1 - 0. 36(K  /K  )3 

[S-'j  a _J_X_  (15) 

1 

0.5  + 0. 14(K  /K  )3 

t4'1!.,  a r-2-  1 06) 

1 

When  K^>  K^,  interchange  all  1 and  2 subscripts  on  both  sides  of 
Equations  (15-16) . 


III.  FAST  MODE  LOCK-IN 

The  preceding  section  assumed  that  the  e and  y motions  had  constant 
phase--that  they  were  "locked-in"--with  the  angle-of-attack  plane.  We 
note  from  Equations  (6-8)  that  the  frequency,  damping  and  roll  effect 


mode  is  much  greater  than  that  of  the  slow  mode.  This  suggests  a much 
simpler  theoretical  model,  namely,  one  in  which  the  e and  y motions  are 
locked-in  with  the  plane  of  the  fast  mode.  Then  Equation  (10)  is  replaced 
by 


0 = <*! 

' 1 

Under  this  assumption,  Equations  (4)  and  (5)  very  quickly  yield 


(17) 


A^1  = “^rS^y/  C0S  *y  + BeV  C°S  [2It' V"LX0] (18) 

(19) 


A$  = 0 
2 


K 

1 


X K + $ [B  Y sin  f + 
11  1 Yi  VY 


BeV  sin  ♦en2ItVLxo] 


(20) 


K = X K 

2 2 2 


(21) 
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Equations  (3)  and  (5)  can  be  combined  and  averaged  to  yield 

*xbpb  + *xc  Pc  ~ \ pb 

P 


- $ K [B  Y sin  6 + B $ e sin  <f>  ] (22) 

1 i Yi  Y e l eJ 


IV.  DISCUSSION 

The  ball  fuze  is  a case  where  the  center  of  mass  motion  can  occur 
but  we  can  assume  y = 0 and  = pc.  Equations  (20-22)  become 

\ = \ K1  + mc  xc  £ G sin  [2  Xt  *1  - lx  pb1_1 


Tx  Pb  = \ Pb  - ♦?  K!  mc  xc  £ e sin 

p 1 1 


In  all  flights,  K remained  less  than  120  milliradians  and  no  spin- 
down  moment  could  be  observed. 

In  Figure  2,  K /K  is  plotted  versus  Mach  number  for  the  20mm 
1 1 

shell  T282E1  with  and  without  the  ball  rotor.  For  Mach  numbers  below  2, 
the  exponential  damping  is  well  determined,  but  above  this  Mach  number, 
considerable  scatter  occurs.  Projectiles  whose  fuzes  did  not  have  the 
ball  rotor  have  damping  rates  that  lie  close  to  -7  per  second  while 
those  with  the  rotor  have  values  that  are  as  much  as  9 per  second 
greater  (i.e.  +2). 

This  phenomenon  for  Mach  numbers  near  2 can  be  ascribed  to  the 
action  of  a locking  spring  which  releases  the  ball  when  the  shell  spin 
rate  is  large  enough.  For  a given  twist  gun  tube  and  air  temperature, 
the  initia1  spin  rate  is  proportional  to  the  launch  Mach  number. 


. mu  iiupu^pnip^pppppp 


The  appropriate  parameters  for  the  T282E1  at  Mach  number  3.3  are 
given  in  Table  1.  For  Mach  numbers  between  2.8  and  3.8,  measured  Kj  falls 

between  .05  and  .12.  The  observed  damping  rate  discrepancies  can  be 
explained  by  Equation  (23)  with  values  of  Ie  sin  <t>£  between  0.09mm  and 

0.28  mm.  These  values  are  possible  for  the  actual  tolerances  observed 
in  this  shell. 

The  actual  spin  histories  of  several  T317's  are  given  in  Reference  4 
and  are  repeated  as  Figure  3.  This  figure  also  gives  the  spin  history  for 
three  T347.  The  T347  shell  has  the  same  external  shape,  mass,  and 
moments  of  inertia  but  no  moveable  internal  components.  In  all  observed 
cases,  the  T317  had  i greater  spin  loss  and  flew  to  a lesser  range.  The 
relative  decrements  between  the  range  of  each  T317  shell  and  the  average 
range  of  the  T347's  is  given  in  the  figure.  Thus  a spin  loss  of  almost 
70  Hz  was  observed  for  a projectile  that  flew  11%  short  of  its  proper  range. 

Unfortunately,  measurements  of  angular  motion  were  not  made.  A 
range  loss  of  11%  would,  however,  require  an  average  angular  motion 
amplitude  of  10°  to  15°.  The  shell  internal  construction  is  fairly 
complicated  but  can  be  theoretically  approximated  by  a single  ring  freely 
sliding  on  a central  shaft  (Figure  4).  The  tolerances  are  quite  small 
but  are  sufficient  to  allow  the  ring  to  cant  at  an  angle  as  large  as 
.004  radian.  When  the  ring  is  fully  canted,  its  center  of  mass  is  on 
the  axis  of  the  shell  and,  hence,  e is  zero.  For  simplicity  we  will 
assume  that  the  ring  is  spinning  with  the  shell.  The  fast  mode  damping 
rate  as  given  by  Equation  (20)  becomes: 

■ \\  * V'xA  - ItcVY  51,1  *Y  [2It+,  - W1  C26) 

The  spin  equation  (Equation  (22))  now  reduces  to 

Vb  = \p  ?b  - Vl(Ixcpb  ‘ 'tcVY  si"  b f27) 

For  a dynamically  stable  shell,  X ^ is  negative.  If  <J>^  is  in  the 

first  or  second  quadrant,  the  second  term  in  Equation  (26)  is  positive 
and,  hence,  K should  grow  to  some  equilibrium  value.  This  equilibrium 

value  of  K can  then  be  used  in  Equation  (27)  to  give  a spin-down  moment. 

Our  theoretical  model  of  the  behavior  of  the  T317  is  based  on  the 
assumption  that  it  has  values  of  A , y and  ^ such  that  its  fast  mode 

motion  grows  to  10-15°  and  causes  the  observed  range  loss  and  spin-down. 


16 


1 


TABLE  1.  PARAMETERS  FOR  T282E1  AT  MACH  NUMBER  3.3 


l = 2.0  cm 


m = 96.6  gm 

I =53.4  gtn/cm2 
x 

I = 394  gm/cm2 


m = 3.0  gm 
c 


CD  . 0.41 

C,  =2.6 

a 


C =2 
LM 

a 


CM  + C = -5.3 
M M* 

q a 


x = 2.9  cm 
c 


V = 1140  m/s 
o 


P = 1900  Hz 
o 


Cu  = 0.15 
M 

pa 


**********  mm 


As  a check  on  this  conjecture.  Equations  (3-4)  were  coded  for  a 
digital  computer.  Although  the  complete  set  of  aerodynamic  coefficients 
are  not  well  known  for  this  shell,  nominal  values  were  used  and  are 
given  in  Table  2.*  Figures  5 and  6 give  the  yawing  motion  and  spin  for 
the  completely  rigid  shell  (y=0) . The  yawing  motion  shows  the  usual 
small  amplitude  slow  mode  limit  cycle  which  is  frequently  observed. 

Computer  runs  were  then  made  for  y = .004,  $ = 45°.  Figure  7 

shows  a rapid  growth  of  the  fast  mode  angular  motion  to  about  18°  and 
a decay  on  the  down  leg  of  the  trajectory.  A range  loss  of  11%  was 
computed  and  a large  spin-down  is  shown  in  Figure  6.  The  computed 
spin-down  is  as  large  as  the  observed  spin-down  but  different  in  detail. 
Thus  we  would  assume  that  the  actual  angular  motion  grew  much  slower  than 
our  computed  motion  but  reached  a much  larger  maximum  value. 

In  summary  then,  the  theory  gives  good  qualitative  agreement  with 
the  observed  behavior  of  the  T317.  In  view  of  the  incomplete  information 
on  its  aerodynamic  properties  and  angular  motion,  better  quantitative 
results  can  not  be  expected. 


V.  CONCLUSIONS 

A theory  has  been  developed  for  the  motion  of  a projectile  with 
a moving  internal  component  that  performs  either  a forced  center  of 
mass  motion  or  a forced  precession  of  its  spin  axis. 

This  theory  gives  a good  explanation  of  the  observed  reduced 
damping  of  shell  with  a ball  rotor  in  its  fuzes  and  the  reduced  range 
and  rapid  spin-down  of  a projectile  with  a ring  on  a central  column. 
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Equations  ( 6-7,  26)  no  longer  apply  hut  quasi-linear  forms  of 
these  equations  aan  be  obtained  by  the  usual  techniques . 6 
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TABLE  2.  PARAMETERS  FOR  T317 


I /T  = .1352 
x t 

I /I  = .0472 
xc  t 


1.34M 
. 125M 


.0251 


J4.4 

Cu  =<1.32  + 3.5M 
Ma  (5.17 


f- .26  + 962 


M<  .88 
. 88<M$1 . 1 
M>1 . 1 


M<  .88 
.88<M«1.1 
M>1 . 1 


<5^0.1 


20mm  Shell  T282E1  with  Arming  Ball  Rotor 


Figure  2.  Fast  Mode  Damping  Rate  for  the  20mm  T282E1 


Simplified  Schematic  of  the  T317 


Computed  Angular  Motion  for  the 


CONOS  ) 
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LIST  OF  SYMBOLS* 


Y 

B . 
YJ 


(1/2) p S£2VC£ 

P 

(1/2) P S^2V[PbCM  + (Lx0/m£2)  C ] 

P«  a 

Cl/'2)p  S£2V[Cm  + CM<  - (ItM2)  CL  ] 
q a a 

(l/2)pS£V2  CM 

a 

Vc  - !tc5 

Vc  - ’tc+j  • J ■ !-2 

m x i 
c c 

(B^.y  cos  ^ + B£i.e  cos  <)>e)[e10].  , j = 1,2 
roll  damping  moment  coefficient.  Equation  (A39) 
lift  force  coefficient,  Equation  (A41) 

Magnus  moment  coefficient.  Equation  (A40) 
damping  moment  coefficients,  Equation  (A40) 
static  moment  coefficient,  Equation  (A40) 


| *In  this  List , the  words  " projectile "body"  and  "internal  component" 

are  used  in  the  sense  that  the  projectile  consists  of  an  external  symmetric 
body  and  an  internal  symmetric  component. 
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LIST  OF  SYMBOLS  (Continued) 


-*■  -;  -*• 
e , e , e 
1 2 3 


unit  vectors  in  the  nonrolling  aeroballistic  system, 
along  the  body's  axis  of  symmetry 


e , e e 
1C’  2C.  3C 


r 

e !i 


unit  vectors  in  the  internal  component's  coordinate 

system,  e along  the  internal  component's  axis  of 

-»• 

symmetry  and  e^c  normal  to  e^,  Equations  (A23-A25) 


i 2ir  i(8-*.)  « 

27  f e J d*  , j - 1,2 


E.(k) 


complete  elliptic  integral  of  the  i-th  kind, 
i = 1,2,  with  modulus  k 


1 » I* 
x’  t 


axial  and  transverse  moments  of  inertia  of  the 
projectile 


*xb’  *tb 


axial  and  transverse  moments  of  inertia  of  the 
body,  Equation  (A20) 


lxc’  *tc 


axial  and  transverse  moments  of  inertia  of  the 
internal  component,  Equation  (A31) 


K. 

J 


length  of  the  j-th  modal  arm,  j * 1,2 


i 

t 


reference  length 


angular  momentum  vector  of  the  projectile 


£ , £ 

b’  c 


angular  momentum  vector  of  the  body  and  internal 
component,  respectively 


xo 


I . p,  + I p 
xb  rb  xc  *c 


m 


the  mass  of  the  projectile,  m^  + 


m 


% 


mass  of  the  body 
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LIST  OF  SYMBOLS  (Continued) 


mass  of  the  internal  component 


components  of  the  aerodynamic  moment  vector  in  the 
nonrolling  aeroballistic  system.  Equations  (A39-A40) 


components  of  the  projectile's  angular  velocity  in 
the  missile-fixed  coordinate  system,  Equation  (Al) 

components  of  the  projectile's  angular  velocity  in 
the  nonrolling  aeroballistic  system,  Equation  (A4) 

the  roll  rate  of  the  body 

(Jc  + 8,  the  roll  rate  of  the  internal  component 

q + ir 

distance  of  the  body  submass  dm  from  the  e -axis, 
Equation  (A12)  1 


vector  from  the  center  of  mass  of  the  projectile  to 
the  center  of  mass  of  the  body,  Equation  (A8) 


vector  from  the  center  of  mass  of  the  body  to  the 
submass  dm  of  the  body,  Equation  (A12) 

distance  of  the  internal  component  submasa  dm 
from  the  e^-axis,  Equation  (A29) 

vector  from  the  center  of  mass  of  the  projectile  to 
the  center  of  mass  of  the  internal  component, 
Equation  (A9) 
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LIST  OF  SYMBOLS  (Continued) 


vector  from  the  center  of  mass  of  the  internal 
component  to  the  submass  dm  of  the  component, 
Equation  (A29) 


reference  area 


S. 

J 


• 10 

(By;j  y sin  + B£^e  sin  4>e)  [e  ].,  j = 1,2 


time 


V 

x,  y,  z 


V V xc 


magnitude  of  the  velocity  vector 

missile-fixed  coordinates,  with  the  x-axis  along 
the  body's  axis  of  symmetry,  the  y-axis  initially 
pointed  down  and  rolling  with  the  body  and  the  z-axis 
determined  by  the  right-hand  rule 


e^- component  of  r^,  and  ?c>  respectively 


e component  of  itc 

angles  of  attack  and  sideslip  in  the  nonrolling 
aeroballistic  system 


the  cant  angle:  the  angle  (assumed  small  and 

constant)  between  the  axis  of  symmetry  of  the 
internal  component  and  the  axis  of  symmetry  of  the 
body 


LIST  OF  SYMBOLS  (Continued) 


— K -K  .)E  k)  + K +K  )E  (k)  , j = 1,2 
ttK?  J 3-J  1 12  2 

3 

j = 1,2 

e.  + e , the  nondimensional  radius  of  Che  circular 
b c 

motion  performed  by  the  center  of  mass  of  the  internal 
component  about  the  axis  of  symmetry  of  the  body 


(~¥)e 


B£  0 e exp  (i8£) 


the  orientation  angle  of  the  angle-of-attack  plane 

i 0 

with  respect  to  the  aeroballistic  axes  (6e  =C); 

in  most  relations,  0 can  be  approximated  by  if) 


the  orientation  angle  of  the  cant  plane  with  respect 
to  the  aeroballistic  axes 


the  polar  angle  of  the  circular  motion  performed 
by  the  center  of  mass  of  the  internal  component 
about  the  axis  of  symmetry  of  the  body 


(A  + Apa)/(2It  - Lx0),  the  j-th  damping  rate 
(Kj/Kj)  for  a rigid  projectile,  j = 1,2 


8 + ia,  the  complex  angle  of  attack 


♦ 


b 


LIST  OF  SYMBOLS  (Continued) 

the  polar  angle  in  the  e -e  plane  of  the  body 

2 3 

submass,  Equation  (A12) 


4>c  the  polar  angle  in  the  e c-e3C  plane  of  the 

internal  component  submass.  Equation  (A29) 


<i>j0  + 'J’jt,  the  orientation  angle  of  the  j-th  modal 
arm,  j = 1,2 

initial  orientation  angle  of  the  j-th  modal  arm, 

j = 1,2 

frequency  of  the  j-th  modal  arm,  j = 1,2  (it  is 

assumed  that  | >$  , that  is,  the  1-arm  is  the  fast 
arm)  1 2 

frequency  of  the  j-th  modal  arm  for  a rigid 
projectile,  j = 1,2 


0-0,  assumed  constant 
Y 

0 -0,  assumed  constant 
e 

♦ -♦ 

2 1 

angular  velocity  of  the  projectile.  Equations  (Al,  A4) 


angular  velocity  of  the  missile-fixed  coordinate 
system,  Equations  (A2,  A4) 

angular  velocity  of  the  nonrolling  aeroballistic 
coordinate  system,  Equations  (A3,  A5) 
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Subscripts 

f 'av 

b,c 

x,y,z 


LIST  OF  SYMBOLS  (Continued) 


body,  internal  component 


vector  components  in  the  nonrolling  aeroballistic 
system 


Special  Notation 

a + ib  complex  representation  of  the  vector  ae  + be 


( ) 


d(  )/dt 
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APPENDIX  A.  DERIVATION  OF  EQUATIONS  OF  MOTION 


A missile-fixed  coordinate  system  can  be  defined  with  x-axis 
along  the  missile  axis  of  symmetry,  z-axis  initially  pointed  down  and 
rolling  with  the  missile,  and  a y-axis  determined  by  the  right-hand 
rule.  If  the  angular  velocity  of  the  missile  in  these  coordinates  is 


<*»  = (p»  q»  r) 


(Al) 


the  angular  velocity  of  the  coordinate  system  is 


(P.  q» 


r) 


(A2) 


For  the  motion  of  symmetric  missiles,  a much  more  convenient  co- 
ordinate system  is  the  aeroballistic  system,  which  pitches  and  yaws 
with  the  missile  but  has  zero  roll  rate.  Its  angular  velocity  vector 
in  missile-fixed  coordinates  is 

ft  = (0,  q,  r)  (A3) 


Transverse  components  of  a vector  in  the  nonrolling  coordinates  will 
be  identified  by  tilde  superscripts.  In  aeroballistic  coordinates, 
the  three  angular  velocity  vectors  assume  the  form 


u = ft  f = p + q e2  + r e3  (A4) 


ft  = q e2  + r e3  (A5) 


where  e , 


are  unit  vectors  along  the  aeroballistic  axes. 


We  now  consider  the  projectile  to  consist  of  its  external 
symmetric  body  with  mass  m^  and  an  internal  symmetric  component  of 

mass  mc  which  is  free  to  move  perpendicular  to  the  body's  axis  of 

symmetry.  The  axis  of  symmetry  of  the  internal  mass  is  assumed  to 
cant  at  constant  angle  y with  respect  to  the  body's  axis  of  symmetry 
and  the  plane  of  this  cant  angle  is  assumed  to  maintain  a constant 
phase  angle  <f>  with  the  angle-of-attack  plane.  If  0 and  6 are  the 
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orientation  angles  of  the  angle-of-attack  plane  and  the  cant  plane, 
respectively,  then 


0 = 0 + <b 

y y 


The  motion  of  the  center  of  mass  of  the  internal  component  will 
be  assumed  to  be  a circular  motion  of  amplitude  £e  and  constant  phase 


angle  with  respect  to  vhe  angle-of-attack  plane.  If  0^  is  the 


phase  angle  of  this  center  of  mass  motion,  then 


0 =0  + 4 

e Te 


With  respect  to  the  body-plus-intemal -component  center  of  mass,  both 
the  body  and  the  component  are  performing  circular  motions  with  radii 
£efe  and  £ec>  respectively.  The  vectors  locating  these  two  centers 

of  mass  with  respect  to  the  projectile  center  of  mass  are 


r,  = x.  e,  + £e,  (e„  cos  0 + e,  sin  0 ) 

b d 1 b v 2 e 3 eJ 


r = x e - £e  (e„  cos  0 + e,  sin  0 ) 

c c 1 c 2 e 3 e' 


x,  m,  + x m =0 
o d c c 


where 


e,  = m e/m 
b c 


ec  = m^  e/m 


m = m,  + m 


The  angular  momentum  vector  of  the  body,  can  now  computed 
from  its  definition8  in  terms  of  a large  number  of  small  submasses, 


8.  Herbert  Goldstein , Classical  Meohanics , Reading , Massachusetts, 
Addison-Wesley  Publishing  Company,  Inc.,  1950. 


if 
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dm,  with  position  vectors  + 


• • 


lb  * |f?b  * V x (?b  * V 


dm 


l«b  rb  X rb  * J fSb  X 


dm 


(All. 


The  rotational  symmetry  of  the  body  can  be  best  exploited  by 
expressing  the  position  vectors  of  the  submasses  in  cylindrical  co- 
ordinates : 


^b  s ?i  + rb  [e2  cos  ^ + e3  sin  (A12; 

The  rotational  symmetry  implies  several  useful  integral  relationships 


i *b  \ cos  *b  (ta 1 1 *b  \ sin  *b da  * 0 


(A13) 


j r£  sin  <j>b  cos  $b  dm  = 0 


(A14) 


j r2  cos2  $b  dm  = J r2  sin2  $b  dm  = (1/2)  j r2  dm  (A15; 


Since  the  body  is  spinning  with  respect  to  the  aeroballistic  axes,  it 
should  be  emphasized  that  $>b  is  not  zero.  The  derivatives  of  rb  and 

l Rb  can,  however,  be  easily  computed  using  the  relation 

i 


i 


i 

I 39 


The  first  cross  product  in  the  angular  momentum  equation 
(Equation  (All))  can  now  be  computed  from  Equations  (A8)  and  (A17)  and 
simplified  by  neglecting  terms  involving  e£  . 


r^  X ?b  = - i eb  § (e2  cos  0£  + e3  sin  0£) 


+ ?bx(^x?b)  ' . I 

4 

-+•*  I 

- ^ i eb  (q  cos  0£  + r sin  0g)  ex  « 


+ (q  x2  - ^ £ eb  0 cos  0£)  e2  ! 

t 

l 

* & *b  " *b  1 Gb  ® sin  V ®3  (A19) 

/ 

/ 

The  integral  in  Equation  (All)  can  be  simplified  quite  easily  to  a ! 

familiar  form.  This  is  < ue  to  the  fact  that  the  aeroballistic  axes 
are  normal  axes  of  inertia  fo-r  the  body. 


where 


‘xb  - J *b  * 


■tb  - / * W'2’  £b> 


Equations  (A19-A20)  can  now  be  used  to  compute  the  angular 
momentum  of  the  body  from  Equation  (All) : 


t = L v £ + L-,  e„  + L~  e 
b xb  1 yb  2 zb  3 


where 


Lxb  ■ Pb  !xb  * “s  % « C0S  \ * f Sin  V 


Lyb  " < "tb  * ”b  Xb)  * Bc  £b  5 C0S 


Lib  * J (!tb  * "b  # * Bx  Eb  B sln  6, 


Pb  ■ *b  1 Be  *c  xc  1 


The  angular  momentum  of  the  internal  component  can  be  computed 
in  a similar  way 


L = m r X r + ( it  X ft  dm 

c c c c I c c 


(A22) 


The  derivative  of  this  vector  is  somewhat  more  complicated  than  the 
derivative  of  R,  since  a term  involving  the  precession  rate,  0,  is 
present.  D 


RC  = ■ Jc  *c  (e2c  sin  *c  - e3C  C0S  ^ 


0 < [v  x sin  0 + r sin  (<f>  + 0 }]e„ 

1 1 ' c y c c y 2 


- [Y  x cos  0 + r cos  ($ 

l i c Y c 


+ 6 )]e,  i 

c y 3 f 


•>  -* 
+ fl  X R 


(A30) 


The  integral  in  Equation  (A22)  can  now  be  computed  in  a similar 
fashion  as  was  done  for  the  integral  in  Equation  (All)  since  symmetry 
conditions  like  Equations  (A13-A15)  apply.  An  additional  term  in- 
voicing the  precession  rate,  §,  must,  of  course,  appear. 


[ R X $ dm  = 4> 

J c c Tc 


I 

xc 


+ 0 [I  e + y (I  -I*.  ) (e„  cos  0 + e,  sin  0 )j 

1 xc  1 xc  tc' v 2 Y 3 Y 

+ f X (n  X itc)  dm  (A31) 
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where 


«*  t 
rc.  1 

5 


* > 


V 


The  integral  on  the  left  side  of  Equation  (A31)  can  easily  be  handled 
if  the  unit  vectors  along  the  normal  axes  of  inertia  of  the  internal 
component  are  used. 


j Rc  x (0  X Ic)  dm  = Qlc  Ixc  elc  ♦ Itc  (a2c  J2c  * n3c  ?3c) 


= (I  - I.  ) n,  e + I S 
v xc  tc'  lc  1c  tc 


(A32) 


where 


Since  the  first  cross-product  in  Equation  (A22)  can  be  handled 
by  an  equation  similar  to  Equation  (A19) , Equations  (A31-A32)  can  be 
used  to  yield  the  angular  momentum  of  the  internal  component  in  the 

aeroballistic  axes.  s 

. * 


t 


= L e, 
xc  1 


+ L-  e + L-  e„ 
yc  2 zc  3 


(A33) 
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where 


;j  l 


L = I p + (I  - I.  ) y (q  cos  8 + r sin  0 ) 

xc  xc  *c  v xc  tc'  n y Y 


+ B e (q  cos  0 + r sin  0 ) 
e c n e e 


Lyc  ' (Itc  * mc  X'P  5 * \ Y cos  \ 


!*i 


+ B 0 e cos  0 
e c e 


L~  * (I.  + m x2)  r + B Y sin  6V 

zc  v tc  c c'  Y Y 


\ ' 


*V  J 


I 


+*4 


+ B 

e 


\ 

0 e_  sin  0_ 

C G 


pc  ‘K + ® 


B 

Y 


b 


Equations  (A21)  and  (A33)  can  be  added  to  yield  the  total  angular 
momentum  of  the  projectile: 


L = 


L e + L~  e + 
x 1 y 2 


L. 

z 


(A34) 
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where 

L = L + (I  - I*  ) Y (q  cos  0 + r sin  6 ) 
x xo  v xc  ter  n y y 


+ b e (q  cos  6 + r sin  0 ) 

e n e e 

L~  = I q + B y cos  0 + B 0 e cos  0 

y t n y Y e £ 

«# 

L~  ■ I.  r + B y sin  0 + B 0 e sin  0 

z t Y Y e e 

ScO  ~ *xb  + Xxc  Pc 
:t  ■ :tb  * rtc  * *b  *b  * mc  xc 

The  differential  equations  for  the  angular  motion  can  now  be 
computed  in  the  usual  way. 

= L e,  + L~  e,  + L~  e,  + ft  X L 

x 1 y 2 z 3 

. = M e,  + M~  e + M-  e„  (A35) 

xl  y 2 z 3 


The  roll  equation  is  obtained  from  the  first  component  of 
Equation  (A35) : 


L + B (q  sin  0 - r cos  0 ) y 

x y n Y Y 


+ B^  0 (q  sin  0£  - r cos  0£)e  = (A36) 
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Since  our  primary  interest  in  the  rolling  motion  is  in  its  average 
behavior,  can  be  replaced  by  its  average. 


(L  ) 
v x'av 


= L 


xo 


- Xxb  P 


, + I p 
b xc  rc 


(A3  7) 


The  transverse  angular  motion  equation  is  obtained  by  multiplying 
the  third  component  of  Equation  (A35)  by  i and  adding  the  result  to 
the  second  component: 


I 


t 


« 

Q 


1 Lxo 


q + r + E = 


M~  + iM~ 
y z 


(A38) 


where 


Q = q + i r 


r = y exp  (i  0^) 


E » B£  0 t exp  (i  0£) 


The  usual  linear  aerodynamic  moment  can  be  written  in  the  form6 


M 


x 


M~  + i M~  = 
y z 


= (1/2)  p S £ V2  C£  (pb  £/V) 
P 


(1/2)  pSIV2 


(p  I/V)  C„  - i CM 
pa  a 


+ CM  (Q£/V)  - i CM  (5  £/V)  y 
q a 


(A39) 


(A40) 
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f 


where 


? = 5 + i a = complex  angle  of  attack 
A good  approximate  relation  between  Q and  l is6 


i Q = I + 


p S V 
2 m 


(A41) 


The  small  lift  term  can  be  neglected  when  we  use  this  relation  in  the 
roll  equation,  but  that  term  does  have  a damping  effect  on  the  complex 
angle  of  attack  equation.  Equations  (A39-A41)  can  now  be  used  to 
write  Equations  (A36)  and  (A38)  in  their  final  forms: 


[xb  pb 


+ I 


p 

xc  rc 


- A. 


Pb  + 


(0 


Bj(0  cos  0 + 


a sin  6 ) 
Y 


) 


av 


. £ 
, 0 (0 


+ B 0 (0  cos  0 + a sin  0 ) 

e e e 


e = 0 


(A42) 


av 


It  l - (A + iLx0)C  - (Ao  + iA  ) l ♦ i (f  + E)  . 0 (A43) 


P« 


\ 

5 

j 


i 

! 


where 


(1/2)  p S *2  V C£ 

P 


i 


A = (1/2)  p S £2  V[CM  + CM  - (It/m  Jt2)  CL  ] 
^ q d a 

Aa  = (1/2)  p S £ V2  CM 

a 

V = a/2)  P s t2  Vtpb  CM  * <LxO/m  l2)  CL  ) 

pa  a 
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Although  Equations  (A42-A43)  were  derived  for  a canted  internal 
component  precessing  with  the  angle-of -attack  plane,  they  oan  be 
applied  to  a canted  component  that  is  fixed  with  respect  to  the  body 
of  the  shell.  For  this  case, 


0 = p,  and  4 


•'  Pc  = ?b 


(A44) 


(A4S) 


Thus  Equations  (A42-A43)  also  are  valid  for  the  motion  of  a shell 
with  mass  asymmetry. 


ArPENDIX  B.  QUASI- LINEAR  SOLUTION  OF  PITCHING  AND  YAWING  MOTION 

The  quasi- linear  approach  assumes  that  the  actual  angular  motion 
can  be  approximated  over  sections  of  the  trajectory  by  the  epicycle 
solution  of  the  linearized  equation  with  constant  coefficients. 


5 = K e 
1 


Hi  + K eH2 

2 


4 . - <ti . t + 6 . 

J J jo 


Equation  (Bl)  can  be  substituted  in  Equation  (A43)  and  the  small  terms 
involving  omitted: 

It(-y|  * 2i  ♦ , V - (Aq  * * i K,  *,) 


- (A  + iA  ) K 

a P i 
*a  L 


= -i(f  + Eje-11^1  + {}  e1^ 


where  $ - -<J>  and  {}  is  the  same  expression  as  on  the  left  side 

2 1 

of  Equation  (B2)  with  the  1 subscript  replaced  by  a 2 subscript. 

This  second  term  on  the  right  side  of  Equation  (B2)  is  periodic 

with  an  average  value  of  zero  while  the  left  side  is  constant.  We 

therefore  have  to  average  the  first  term  on  the  right  side  to  obtain 

its  contribution  to  the  quasi-linear  frequencies  and  damping  rates. 

9 and  0 in  F and  E respectively  can  be  computed  from  0 using  the 
Y £ 

definition 


[K  eX<t>1  + K e1<i>2]  [K2  + K2  + 2K  K cos 
1 2 12  12 


When  exp  (i6)  is  multiplied  by  exp  (-i<J>^),  we  see  that  the  product 
is  a periodic  function  of 


Now 


2tr 


ICe 


ie 


= ■< 


li  eie  e-^l] 


2it 


2tt 


+ 1 


i*/  e1t0"<l>l)d$l/2iT 


•1  r i^i 

i4»1  [e  ] 


r 


[ C9ei0) * e_i^ 1 ] av  =-|  [ieei6e”i^1] 


2rr 


2tt 


| t (e*0)  e ^Jd$  r/2ir 
vo 


•12  r 

[e  ] 


where 


2jt 


[e 


? 3 d4> 


(B4) 
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Equations  (B4-B5)  can  be  used  in  averaging  Equation  (B2)  and  the 
result  separated  into  real  and  imaginary  parts: 


r \2  - (L  - C K_1H  + A = 0 
t 1 XO  1 1 a 


(B6) 
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(B7) 


t' 


l f ; 

I i* 

u\ 


K = X K + l S [2lJ  - L 1 

1 1111  * 1 xO 


-1 


where 


C.  = [B  .y  cos  d>  + B |.e  cos  $ ] [e10] . 
3 YJ  Y e ] e J 


S.  = [B  .y  sin  <{>  + B 6.e  sin  <p  ] [e*0] . 
J yj  Y e 3 TeJl  J3 


P 

W'. 

ill 


l r 


I k,i 

i i 

h 

V , 


i 


t1 

to 

»K 


B . = I p - I.  $ . 
Y3  xcrc  tcT3 


X.  = [A  J.  + A ] [21. - L J 

] 1 qy3  paJ 1 tj  xOJ 


-1 


Ip  a similar  fashion,  the  frequencies  and  damping  rates  for  the 
second  mode  can  be  derived. 


• 1. 


*22  - <LX0  - C2  K2  > +2  * A«  ■ ° 


K =11  * i S [21  ♦ - L .] 
2 22  22  *2  x0 


-1 


Let  A4>j  * ^ where  is  the  j-th  frequency  of  the  rigid 

body  consisting  of  the  external  body  and  the  internal  component. 
Equations  (B6)  and  (B8)  become 


L ($?  + 2 l.  A$.)  - L . (i,r  + A *.) 

t VTjr  jr  3 xo  3r  1 


+ C.  K.'1  + A = 0 

J J jr  a 


fl 


Thus 


A$ . = -i-J-  Vr- 


J 


2It‘*>jr  " LX0 


1 «.•? 
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(B8) 

(B9) 

(BIO) 

I 

CBll) 
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APPENDIX  C.  QUASI- LINEAR  SOLUTION  OF  THE  ROLL  EQUATION 

The  average  of  the  periodic  y-term  in  the  roll  equation  (Equation 
(A42))  can  be  obtained  from  the  following  relations*: 


[sin  (0  70-)] 

1 v y j * 1 av 


i[(e-<M+$  ] ] 

e J Y -e  J ' 


, i0-,  , 

= [e  1 . sin  4 
L Y 


ii  i -i0 

[6  cos  0 + a sin  0 ] = R{£  e 

1 y y av  av 


= UK  sin  (0-4)  + 4 K sin  (0  -4  )LV 
11  Y 1 2 2 Y 2 av 


C4>  K [e10]  + 4 K [e10]  ) sin  4V 

1 1 1 2 2 2 Y 


[0(0  cos  0 + a sin  0 )]  = [(0  sin  0 - a cos  0 ) 

i w y yJ1av  Y Y 


■ 0 sin  0y  + a cos  0y]av 

■j  -ie 

i«  • Y).„ 


[42  K sin  (0  -4  ) + 42  K sin  (0  -4  )] 

1 1 1 Y V 2 2 Y 2 Jav 


(42  K [e10]  + 42  K [e10]  )sin  4 

1 1 1 2 2 2 Y 


Equations  (C2-C3)  and  similar  relations  for  subscript  e can  now  be 
used  to  average  Equation  (A42))  which  reduces  to 


I , P,  + I P = A„  P,  - 4 K S - 4 K S 
xb  b xc  c £p  b i i i 2 2 2 

*For  the  Q's  of  this  report , f (,0"!*1l/)  c 


Note  that  the  spin  moment  depends  on  y sin  <f>  and  e sin  4>  , the 

i$Y  i^e  ^ e 

out-of-plane  components  of  ye  T and  ee  . Moreover,  only  the  out- 
of -plane  components  appear  in  the  damping  rate  equations  (B7)  and  (B9) . 
If  S and  S2  are  positive,  their  effect  is  to  undamp  the  fast  (K^  mode. 


damp  the  slow 


(K2)  mode  and  reduce  the  spin  rate. 
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